Celestial Mechanics and Dynamical Astronomy, Vol.94. No.l. (2006) pp. 1-15. 

Gravitational Field of 
Fractal Distribution of Particles 

Vasily E. Tarasov 

Q ' Skobeltsyn Institute of Nuclear Physics, 

(N 

i-i ' E-mail: tarasov@tlieory.sinp.msu.ru 

< 

CN ! Abstract 

In this paper we consider the gravitational field of fractal distribution of parti- 
cles. To describe fractal distribution, we use the fractional integrals. The fractional 
—1. I integrals are considered as approximations of integrals on fractals. Using the frac- 

'^ ' tional generalization of the Gauss's law, we consider the simple examples of the 
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■ ■■ • 1 Introduction 

The aim of this paper is to consider the general properties of gravitational field that is 
generated by a fractal distribution of particles. This problem is nowadays particularly 
relevant. In fact, there is a general agreement that galaxy distribution exhibits fractal 
behavior up to a certain scale [H ^. The eventual presence of a transition scale towards 
homogeneity and the exact value of the fractal dimension are still matters of debate. 
Moreover it has been observed that cold gas clouds of the interstellar medium has a 
fractal structure, with 1.5 < D < 2 in a large range of length scales [21 E]. Derivatives 
and integrals of fractional order fS| have found many applications in recent studies of 
fractal structures. Fractional analysis can have numerous applications: kinetic theories 
inmiEllII]; statistical mechanics 1311011121; dynamics in complex media [THllTHIT^ITEl lTTj: 



and many others. The new type of problem has increased rapidly in areas in which the 
fractal features of a process or the medium impose the necessity of using non-traditional 
tools in "regular" smooth physical equations. In order to use fractional derivatives and 
fractional integrals for fractal distribution, we must use some continuous medium model 

mmniini. 

We propose to describe the fractal distribution by a "fractional" continuous medium 
[T^ . where all characteristics and fields are defined everywhere in the volume but they 
follow some generalized equations which are derived by using fractional integrals. In many 
problems the real fractal structure of medium can be disregarded and the fractal distribu- 
tion can be replaced by some "fractional" continuous mathematical model. Smoothing of 
microscopic characteristics over the physically infinitesimal volume transforms the initial 
fractal distribution into "fractional" continuous model P3j that uses the fractional inte- 
grals. The order of fractional integral is equal to the fractal dimension of distribution. 
The fractional integrals allow us to take into account the fractality of the distribution. 
Fractional integrals are considered as an approximation of integrals on fractals |18| IT^ . 
It was proved in |TH] that integrals on net of fractals can be approximated by fractional 
integrals. In Ref. [H], we proved that fractional integrals can be considered as integrals 
over the space with fractional dimension up to numerical factor. In order to prove, we 
use the formulas of dimensional regularizations [201 • 

In this paper, we consider gravitational field of fractal distribution of particles. Frac- 
tal distribution is described by fractional continuous medium model [HI Uni El EI- In 
the general case, the fractal distribution of particles cannot be considered as continuous 
medium. There are points and domains that are not filled of particles. In Ref. j2], 
we suggest to consider the fractal distributions as special ("fractional") continuous me- 
dia. We use the procedure of replacement of the distribution with fractal mass dimension 
by some continuous model that uses fractional integrals. This procedure is a fractional 
generalization of Christensen approach [2Ij. Suggested procedure leads to the fractional 
integration and differentiation to describe fractal distribution. 

In Section 2, the density of fractal distribution of mass is considered. In Section 3, we 
consider the simple examples of the gravitational field of homogeneous fractal distribution. 
In Section 4, the examples of gravitational quadrupole moments for fractal distribution 
are considered. Finally, a short conclusion is given in Section 5. 



2 Mass and Balance of Mass for Fractal Distribution 

2.1 Mass of Fractal Distribution 

Let us consider a fractal distribution of particles. For example, we can assume that 
particles with a constant density are distributed over the fractal. In this case, the number 
of particles A^ enclosed in a volume of characteristic size R satisfies the scaling law 

N{R) ~ i?^, (1) 

whereas for a regular n-dimensional Euclidean object we have N{R) ~ i?". 

For distribution of particles with number density n(r, t), we have that the mass density 
can be defined by 

p(r,t) =mn(r,t), (2) 

where m is the mass of a particle. The total mass of region W is then given by the integral 



MiW)= / pir,t)dVs, (3) 

Jw 

or M{W) = mN{W), where N{W) is a number of particles in the region W. The 

fractional generalization of this equation can be written in the following form: 

M{W)= f pD{r,t)dVD, (4) 

Jw 

where D is a mass fractal dimension of the distribution, and dVn is an element of D- 
dimensional volume such that 

dVD = C3{D,r)dVs. (5) 

For the Riesz definition of the fractional integral j3], the function C3{D,r) is defined by 

the relation 

23-^r(3/2) _3 

^^^^'")- TiD/2) '"' • (^) 

The initial points of the fractional integral are set to zero. The numerical factor in Eq. 
© has this form in order to derive usual integral in the limit D ^ {3 — 0). Note that the 
usual numerical factor 'j^^(D) = T{l/2)/2^n^^'^T(D/2), which is used in Ref. j3] leads to 
73-1(3 - 0) = l/[47r3/2] in the limit D ^ {3 - 0). 

For the Riemann-Liouville fractional integral |S^, the function C3{D,r) is defined by 

I T'ilT'l ' 
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Here we use Cartesian's coordinates x, y, and z. In order to have the usual dimensions of 
the physical values, we can use vector r, and coordinates x, y, z as dimensionless values. 
Therefore the density po has the dimension of mass. 

Note that the interpretation of fractional integration is connected with fractional di- 
mension jHl nni ■ This interpretation follows from the well-known formulas for dimensional 
regularizations [201 ■ The fractional integral can be considered as an integral in the frac- 
tional dimension space up to the numerical factor T(D/2)/{27r^^'^T{D)). 

If we consider the ball region W = {r : |r| < R}, and the spherically symmetric 
distribution of particles (n£)(r,t) = n£){r)), then we have 

N{R) = 47r )' ' / nD{r)r''-^dr, M{R) = mN{R). (8) 

i {i-'l 2) Jo 

For the homogeneous {n£){r) = uq) fractal distribution, we get 

/VW-4™„H!^™^^iJ-. (9) 

Fractal distribution of particles is called a homogeneous fractal distribution if the power 
law N{R) ~ R^ does not depend on the translation of the region. The homogeneity 
property of the distribution can be formulated in the following form: for all regions W 
and W such that the volumes are equal 1^(1^) = V{W'), we have that the number of 
particles in these regions are equal too, N(W) = N(W'). Note that the wide class of 
the fractal media satisfies the homogeneous property. In many cases, we can consider 
the porous media [SHI 121], polymers J2SI, colloid aggregates |^, and aerogels ^^ as 
homogeneous fractal media. In Refs. ^1 |^, the continuous medium model for the 
fractal distribution was suggested. Note that the fractality and homogeneity properties 
can be realized for the fractional continuous model in the following forms: 

(1) Homogeneity: the local number density of homogeneous fractal distribution is trans- 
lation invariant value that has the form n{r) = tiq = const. 

(2) Fractality: the number of particles in the ball region W obeys a power law relation 
Nd{W) ~ R^ 1 where D < 3, and R is the radius of the ball. 

2.2 Flow of Fractal Medium 

For distribution of particles with number density ?7,(r, t) flowing with velocity u = u(r, t), 
the resulting density J(r, t) is given by 

J(r, t) = mn(r, t)u, (10) 



where m is the mass of a particle. Measuring the field J(r, t) passing through a surface 
S = dW gives the flow (flux of mass) 



I{S) = ^jiS)= [i3,dS,), (11) 

Js 



where J = J(r, t) is the flow fleld vector, dS2 = dS2n is a differential unit of area pointing 
perpendicular to the surface S, and the vector n = rijtejt is a vector of normal. The 
fractional generalization of this equation for the fractal distribution can be written in the 
following form 



I{S)= /(J(r,t),dS,), (12) 

Js 

where we use 

dS, = C2id,v)dS2, C2id,r) = ,^^\r\''~\ (13) 

Note that 6*2(2, r) = 1 for d = 2. In general, the medium on the boundary dW has the 
dimension d. In the general case, the dimension d is not equal to 2 and is not equal to 
(D-l). 

2.3 Equation of Continuity for Fractal Distribution 

The change of mass inside a region W bounded by the surface S = dW is always equal to 
the flux of mass through this surface. This is known as the law of mass conservation or 
the equation of balance of mass |16j. If we denote by J(r, t) the flow density, then mass 
conservation is written 

™ - -'iS), (14) 

or, in the form 

^ f pD{r,t)dVn = - / (J(r,t),ciS,). (15) 

In particular, when the surface S = dW is flxed, we can write 

Using the fractional generalization of the mathematical Gauss's theorem (see Appendix), 

we have 

/ {3{r,t),dS,)= [ C^\D,r)^(C2{d,r)Mr,t))dVD. (17) 

Jaw Jw (^^k ^ ' 

Substituting the right hand sides of Eqs. (fT^ and (fTTj) in Eq. (I15|) . we flnd the equation 
of balance of mass in differential form 

C3(Ar)^^^^ + ^(C2(rf,r)J,(r,t)) =0. (18) 

5 



This equation can be considered as a continuity equation for fractal distribution of parti- 
cles [in]. 

3 Gravitational Field of Fractal Distribution 

3.1 Gravitational Field 

For a point mass M at position r' the gravitational field F at a point r is defined by 

F = GM ^^^, (19) 

where G is a gravitational constant. 

For a continuous distribution p(r') of mass, the gravitational field F at a point r is 
given by 

F3(r) = G I T^^r^A^'Wi- (20) 



w 



r — r 



The fractional generalization of this equation for a fractal distribution of mass is given by 
the equation 



r' — r 



FdW = G / -,PD{r')dV'r,. (21) 



w 



r — r 



where dV^ = G3{D,r')dV^. Eq. ()2H1 can be considered as Newton's law written for a 
fractal distribution of particles. 

Measuring the gravitational field passing through a surface S = dW gives the gravi- 
tational filed flux 



^FiS)= /(F,rfS2), (22) 

Js 

where F is the gravitational field vector, and dS2 is a differential unit of area pointing 
perpendicular to the surface S. 

3.2 Gauss's Law for Fractal Distribution 

The Gauss's law tells us that the total flux ^f{S) of the gravitational field F through a 
closed surface S = dW is proportional to the total mass M{W) inside the surface: 

^p{dW) = 47tGM{W). (23) 

For the fractal distribution. Gauss's law states 

f{FD,dS2) = 47rG [ po^Vn, (24) 

Js Jw 



where F = F(r) is the gravitational field, and PD(r) is the mass density, dVo = C^^D, r)dV3, 
and G is the gravitational constant. 

Gauss's law by itself can be used to find the gravitational field of a point mass at rest, 
and the principle of superposition can then be used to find the gravitational field of an 
arbitrary fractal mass distribution. 

If we consider the spherically symmetric fractal distribution PD(r) = Pni^), and the 
ball region W = {r : \r\ < R}, then we have 

M{W)=47i pD{r)C3{D,rydr, (25) 

Jo 

where C3{D,r) is defined in Eq. ©, i.e., 

M{W) = 4vr ^^^J^^ > j^ po{r)r^-'dr. (26) 

Using the sphere S' = {r: |r| = _R}asa surface S = dW, we get 

^p{dW) =A7r R^Fd{R). (27) 

Substituting Eqs. (PB|) and (|77|) in the Gauss's law PH|). we get the equation for gravita- 
tional field. As the result, the Gauss's law for fractal distribution with spherical symmetry 
leads us to the equation for gravitational field 

MR) = ''^Tr.Zti^^ r po{r)r^''dr. (28) 



R'T{D/2) J, 

For example, the gravitational field of homogeneous (pD(r) = po) spherically symmetric 
fractal distribution of mass is defined by 

f.W^P. '%™ i^-^fi-. (29) 

4 Multipole Moments for Fractal Distribution 

4.1 Multipole Expansion 

A multipole expansion is a series expansion of the effect produced by a given system in 
terms of an expansion parameter which becomes small as the distance away from the 
system increases. Therefore, the leading one of the terms in a multipole expansion are 
generally the strongest. The first-order behavior of the system at large distances can 
therefore be obtained from the first terms of this series, which is generally much easier 
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to compute than the general solution. Multipole expansions are most commonly used in 
problems involving the gravitational field of mass aggregations, the gravity and magnetic 
fields of mass and fiow distributions, and the propagation of electromagnetic waves. 

To compute one particular case of a multipole expansion, let R = X^efe be the vector 
from a fixed reference point to the observation point; r = Xke^ be the vector from the 
reference point to a point in the distribution; and d = R — r be the vector from a point 
in the distribution to the observation point. The law of cosines then yields 

2 

d' = R'[l + ^^-2^cos9^, (30) 

where d = |d|, and cos 9 = (r, R)/(ri?), so 



/yi ^ /yi 



d = R\ll + ^-2-cos 6. (31) 



Now define e = r/R, and x = cos 6, then 



- = -(l-2« + .=) , (32) 



-1/2 



But ( 1 — 2ex + e^ ) is the generating function for Legendre polynomials Pn{x) as 

follows: 



oo 



l-2ex + e'^ '^' = ^e"P„(x), (33) 

n=0 

SO, we have the equation 



n=0 

The gravitational potential U (F = — Vf/) obeys (l/d) law. Therefore, this potential can 
be expressed as a multipole expansion 

oo „ 

^ = -^J2l^^ r"P„(cos 9)pD{r)dVD. (35) 

n=0 ^ -^W 

The n = term of this expansion can be pulled out by noting that Po{x) = 1, so 



R 

The nth term 



U = -% I PD{r)dVD - G V -^ / r"P„(cos 9)pDir)dVD. (36) 



^^ = -iMn I ^"^n(cos e)pD{r)dVD (37) 

is commonly named multipole, and for n = 2 - quadrupole. 



4.2 Gravitational Moment of Fractal Distribution 

Gravitational moments describe the nonuniform distribution of mass. The gravitational 
quadrupole term is given by 



-K Jw 



--|- f r2(3cos2 e - l)pD{v)dVD = (39) 



The quadrupole is the third term in a gravitational multipole expansion, and can be 
defined by 

k,l=l 

where G is the gravitational constant, R is the distance from the fractal distribution of 
mass, and M^i is the gravitational quadrupole moment, which is a tensor. 
The gravitational quadrupole moment is defined by the equation 



Mki= / i3xkXi-r'5ki)pDir)dVD, (42) 

Jw 

where Xk = x, y, or z. From this definition, it follows that 

3 

Mki = Mik, and J] M^^ = 0. (43) 

fc=i 

Therefore, we have M^z = —M^x — Myy In order to compute the values 



^i?^= [ [2x' - y' - z']po{r)dVn, 
Jw 



(44) 



M!,y^ = [ [-^' + 2|/' - z')]pnir)dVn, (45) 

Jw 

Mif = [ i-x' -y' + 2z')]pD{r)dVD, (46) 

Jw 

we consider the following expression 

M{a, f3,-f)= [ [ax^ + f3y' + -fz')]pD{r)dVn, (47) 

Jw 

9 



where we use the Riemann-Liouville fractional integral ^, and the function C3{D,r) in 
the form 

Cs{D,r) = l^^, a = D/3. (48) 

Using Eq. (jUj), we have 

M(f =M(2,-1,-1), Mif) = M(-l,2,-l), M(f) = M(-l,-l,2). (49) 

4.3 Quadrupole Moment of Fractal Parallelepiped 

Let us consider the example of gravitational quadrupole moment for the homogeneous 
(pD(r) = po) fractal distribution in the parallelepiped region 

W = {{x;y;z) : 0<x<A, 0<y<B, 0<z<C}. (50) 

If we consider the region W in the form ()5()|1 , then we get ()47|1 in the form 

^^«' ^' ^) = (Z^tf2%, (^^' + f^B' + lC% (51) 

The total mass of this region W is 

M(iy)=Po/rfVz. = 4^- (52) 

Jiy a i [a) 

Therefore, we have the following equation 

Mia,f3,^) = ^M{W){aA' + I3B' + ^C'), (53) 

where a = D/3. If Z^ = 3, then we have 

M(a, P, 7) = ^M{W){aA^ + (3B^ + 7^'). (54) 

As the result, we get gravitational quadrupole moments M^^ of fractal distribution in 
the region W: 

MS' = 7^ M»', (55) 

r{3) 



where M^^ are moments for the usual homogeneous distribution {D = 3). By analogy 

if 



with these equations, we can derive M^^ for the case k ^ I. These quadrupole moments 



are 



Using 2 < D < 3, we get the relations 



^If^ = 77^^ ^if ' (MO- (56) 



3D AD"^ 



Quadrupole moment of fractal ellipsoid is considered in Appendix. 
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5 Conclusion 

The fractional continuous models for fractal distribution of particles can have a wide ap- 
plication. This is due in part to the relatively small numbers of parameters that define 
a random fractal distribution of great complexity and rich structure. In many cases, the 
real fractal structure of matter can be disregarded and the distribution of particles can be 
replaced by some fractional continuous model ^^ E] . In order to describe the distribu- 
tion with non-integer dimension, we must use the fractional calculus. Smoothing of the 
microscopic characteristics over the physically infinitesimal volume transforms the initial 
fractal distribution into fractional continuous model that uses the fractional integrals. 
The order of fractional integral is equal to the fractal dimension of the distribution. The 
fractional continuous model for the fractal distribution allows us to describe dynamics of 
a wide class of fractal media [1^1 IT71 I5]. 

The suggested results can have a wide application to galactic dynamics and cosmology. 
In particular, there is strong evidence that the distribution of mass beyond the scale of 
clusters of galaxies is fractal, with D ^ 1.2, corresponding to a power-law two-point 
correlation function with exponent equal to —1.8. Fractal distributions may also be 
present within gravitational systems of a smaller scale, for example, galaxies. However, 
the results are incomplete in the following aspect: it is known that the fractal distribution 
of mass in the Universe is characterized by large density fluctuations, even within the 
fractal volume, therefore it is far from homogeneous. The nonhomogeneity of the fractal 
distribution can be described by the suggested fractional continuous model |14.|. The 
fluctuation deviation from homogeneity can be parametrized by the so-called n-point 
correlation functions, with n > 2 [2HI 123 EOl E] • Such density fluctuations are important 
in that they produce terms in the force field which can be described only statistically. 
An elementary example is the random distribution of particles in a three-dimensional 
sphere. Although the distribution can be considered as uniform when viewed at the scale 
of the sphere (with D = 3), the Poisson noise of the density field will cause gravitational 
clustering in small scales that will finally prevail the overall evolution of the system, i.e., 
the latter will be quite a different from the evolution of a perfectly homogeneous sphere. 
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Appendix 

Fractional Gauss's theorem 

In order to realize the representation, we derive the fractional generalization of the Gauss's 
theorem 

/ {3{r,t),dS2)= f diviJir,t))dVs, (58) 

JdW Jw 

where the vector J(r,t) = Jk^k is a field, and 

^^^(J) = ^ = 1^- (59) 

or oxk 

Here and later we mean the sum on the repeated index k from 1 to 3. Using Eq. (J13p . we 
get 

(J(r,t),rfS,)= / C2(rf,r)(J(r,t),dS2). (60) 



Note that we have 6*2(2, r) = 1 for the d = 2. Using the usual Gauss's theorem ()58|). we 
get 

/ C2id,r){J{r,t),dS2)= [ dtviC2id,r)Jir,t))dVs. (61) 

JdW Jw 

Equations (0) and (0) in the form dV^ = C^^{D, r)dVD allows us to derive the fractional 
generalization of the Gauss's theorem: 

/ {J{r,t),dSd)= [ C,\D,r)dtv(C2{d,r)3{r,t)) dVn. (62) 

JdW Jw ^ ' 

Quadrupole Moment of Fractal Ellipsoid 

Let us consider the example of gravitational quadrupole moment for the homogeneous 
(pD(r) = po) fractal distribution in the ellipsoid region W: 

^+$+i;<i. (63) 

If we consider the region W in the form ()63|) , then we get ()47|) in the form 

^(«'^'^) = it'^t^X ^ i^^'Kii^) + PB'K2{a) + ^C'K,{a)), (64) 

[6a + 2)L ■^[a) 

where a = D/3, and Ki{a) (i=l,2,3) are defined by 

Ki{a) = L(a + 1, a - 1, 27r)L(a - 1, 2a + 1, vr), (65) 
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K2{a) = L(a - 1, a + 1, 27i)L{a - 1, 2a + 1, n), (66) 

K:i{a) = L(a - 1, a - 1, 27v)L{a + 1, 2a - 1, vr). (67) 

Here we use the following function 

,x 2/ r/\ , , ,,„, , ,,^ / r(n/2 + l/2)r(m/2 + l/2) , , 

TT Jo TT l(n/2 + m/2 + lj 

If /} = 3, we obtain 

M{a, P, 7) = |^^^(aA2 + PB' + ^C% (69) 

where we use Ki = K2 = K3 = An/S. The total mass of this region W is 

^ ^ ^ Jiy 3ar3(a) r(3a/2) ^ ^ 

If /} = 3, we have the total mass 

M{W) = po / dVs = —poABC. (71) 

Using Eqs. (jH!I|l and (ffn|) . we get the quadrupole moments (pn|) for fractal ellipsoid 

M(a, /?, 7) = ^^("^' + PB' + lC% (72) 



where a = D/3. If i^ = 3, then we have the well-known relation 

M{a,(3,'y) = ^^^^{aA^ + (3B^ + ^C^). (73) 
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